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Abstract 

The decay process H — > 77 is examined in a model-independent way 
within the effective Lagrangian approach. Contribution of a set of irre- 
ducible one-loop diagrams involving SU(2) x U(\) invariant dimension-six 
bosonic operators is evaluated explicitly. The calculation is intended to 
fill a gap that exists in the current literature on the subject. 
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1 Introduction 



The standard model (SM) of fundamental interactions has so far proved to be 
a very successful theory, at least from the point of view of present-day phe- 
nomenology. The most prominent unsolved issue is the mechanism of the elec- 
troweak symmetry breaking (EWSB) and a possible existence of the related 
Higgs scalar sector. Current precision tests of the SM seem to point towards a 
relatively light Higgs boson, i.e. such that run < 200 GeV or so (see e.g. pQ). 
Needless to say, such an "evidence" is only indirect, as it is based on a success- 
ful fit of calculated effects of the relevant closed-loop Feynman graphs to the 
precision electroweak data. Note also that there is another message from the 
EWSB sector, which is more general and rather spectacular. For the famous 
parameter p — /m 2 z cos 2 9w (with mw and mz denoting the vector boson 
masses and Oyy being the weak mixing angle) the experimental data show that, 
with great accuracy, p = 1. Such a value is naturally obtained from the Higgs 
sector with doublet structure, or, more generally, from an EWSB mechanism 
respecting the custodial SU(2) symmetry 

On the other hand, a general opinion prevailing nowadays is that SM cannot 
be, for various reasons, an ultimate model of particle physics - in other words, 
it should be understood as an effective theory valid below the energy scale of 
the order O(lTeV) and somewhere above that scale, some kind of new physics 
is to be expected (for a review, see e.g. 0). A general framework for describing 
the effects of the physics beyond SM is an appropriate effective Lagrangian. 
This incorporates, besides the standard renormalizable interaction terms with 
dimension four and three, also non-renormalizable higher-dimensional terms 
involving negative powers of a A, the relevant scale of new physics. Thus, 
the generic form of such an effective electroweak Lagrangian reads 

N n (n) 
n>5 i— 1 

When constructing such an extension of the SM, gauge invariance under SU (2) x 
U(l) is required and, in addition to that, one can make a specific assumption 
concerning the A. We shall assume that A ^ v, where v = 246 GeV is the usual 
electroweak scale. This corresponds to the so-called "decoupling scenario" , in 
which the SU(2) x U(l) symmetry is realized linearly (concerning the termi- 
nology, see e.g. 0], [S] and references therein). In particular, this means that 
the unphysical Goldstone bosons enter via the standard complex SU(2) dou- 
blet, along with the elementary Higgs boson field and, generally, the low-energy 
spectrum is supposed to be identical with that of SM. Explicit representation of 
the non-renormalizable higher-order terms invariant under SU(2) x U(l) gauge 
symmetry can be found in the literature, see in particular j^j, [7j. Contributions 
of the higher-dimensional terms to the amplitudes of low-energy processes are 
then in fact suppressed by powers of the ratio v j A. 

In future experimental searches for the Higgs scalars, the process H — > 77 
might play an important role, since a light scalar boson (i.e. that with ran < 
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2m w) could be detected through this decay mode at LHC (cf. [8] and refer- 
ences therein). Of course, there is no tree- level coupling ^77 within SM: the 
lowest-dimensional interaction term of this type is necessarily proportional to 
HF^yF^" (where F^ denotes the electromagnetic field strength) if one is to 
maintain the electromagnetic gauge invariance. However, such a monomial has 
obviously dimension five and is therefore non-renormalizable. Consequently, SM 
yields a calculable (i.e. ultraviolet-finite) result for H — > 77 matrix element at 
one-loop level and this was obtained long ago [Sj • It can also be sensitive to pos- 
sible effects of the new physics described schematically by the expansion and 
in current literature there are several papers dealing with this issue (see e.g. |1(J) . 

[S]); let us remark that contributions of higher-dimensional effective inter- 
actions to the inverse process 77 — * H have also been considered previously [T2"]. 
Nevertheless, what seems to be still missing, is a one-loop calculation involving 
a complete set of SU(2) x U(l) invariant operators that would yield the leading 
correction (in powers of 1/A) to the SM result. Some important steps toward 
this goal have already been made previously in the papers JU|, and [S]. 
In the present paper we would like to fill a gap that, in our opinion, occurs in 
the existing treatments. In particular, our effort is concentrated on the part of 
the calculation that is perhaps most tedious from the technical point of view: 
we evaluate one-particle irreducible purely bosonic one-loop graphs involving 
just one insertion of a dimension-six effective interaction vertex from l|T[l. with 
other ingredients descending from SM. An important consistency check of our 
computation is the transverse tensor structure of the matrix element - this is 
achieved, similarly as in the SM case, only when a set of relevant graphs is 
summed. Thus, the results presented in this paper are of rather technical na- 
ture, but we believe that they may prove useful in a future complete calculation 
based on the SU(2) x [7(1) effective Lagrangian of the type Q. 

The paper is organized as follows. In the next section we give a brief sum- 
mary of the basic SM results for H — > 77. In Section 3, the set of SU(2) x U(l) 
bosonic operators with dimension six is shown explicitly and in Section 4 we 
display the main features of the whole calculation. The results are described 
in some detail in Section 5 and Section 6 contains their short summary and 
some concluding remarks. To make the paper self-contained, we have added 
the Appendix containing the definitions and explicit analytic formulae for the 
Passarino-Veltman functions that are necessary for the evaluation of loop inte- 
grals. 
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2 Standard Model Results 



To begin with, let us reproduce here - for reference purposes - the form of the 
electroweak SM Lagrangian. In the £/-gauge it reads (for a concise summary 
see e.g. fTS] 1 ) 

^int WS) = E Q f^h^K + J?CC + ^NC 
f 

+ ig(W°W-d t "W + " + W-W+d^W^ + W+W^d^W-") 

- g 2 [\{W- ■ W+f - \{W-f{W + f + {W°f(W- ■ W + ) 

- (W~ ■ W°)(W+ ■ W )] 

+ 9 m w W-W^H + ^- 9 m z Z^H 

- E l^ffH - - ^<H\ (2) 
t-T 1 2 m w 4 m w 32 

where we have denoted 

VK° = cobBwZh + bwBwAh. (3) 

The Jifcc an d ^£nc are the usual symbols for charged and neutral current 
interactions of fermions with vector fields; an explicit form of these terms will 
not be needed for our calculations. 

As we have already stressed in the Introduction, for obvious dimensional 
reasons there is no direct tree-level .H77 coupling within the SM Lagrangian. 
The one-loop SM graphs for the decay process H — > 77 are shown in FigQ] 
Note that at this level, the fermionic part is dominated by the contribution of 




Figure 1: SM one- loop diagrams for H — > 77. Crossing of the external photon 
lines is implicit in contributions of the graphs (a) and (c). 
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the heaviest fermion species. Therefore we consider here the top quark loop 
only. 

A straightforward calculation yields the following result for the decay matrix 
element: 

M = (F VB + F top )[g^ - —P-y a (py^k), (4) 

where the "formfactors" -Fvb and -Ftop represent the contributions of the vector 
boson and top quark loop respectively. These are expressed through 

FyB = 8?r 2 ^ + 6rr?~ + ( 2m ^ _ m H) c o(m 2 H , 0, 0, m^, m^, m%/)\ , (5) 

2 

Ftop = ^^[2 + (4m 2 - m 2 H )Co(m 2 H ,0,0,mlm 2 t ,m 2 t )] , (6) 

with Co denoting the relevant Passarino-Veltman function (see e.g. H3|); the 
definitions and explicit analytic formulae for a set of these functions are sum- 
marized in the Appendix. 

Notice also that in writing the result we have pinpointed the tensorial 
factor 

TP .v = gtl v _ 2m - 2 k> 1 p v , (7) 

which is transverse with respect to the external photon momenta, i.e. 

PfjT^" = 0, KT^ = (8) 

as anticipated on the basis of the electromagnetic gauge invariance. The appear- 
ance of the transverse tensor structure Q will serve as an important consistency 
check also in our subsequent calculations with the effective Lagrangian Q). 

Finally, let us add that from the matrix element (@J one obtains the decay 
width 

. . 8na 2 GF fnlv , „ ,-, ,9 

r SM (i?^77) = ^^iFtop + .FVBl 2 , (9) 

(cf. m). 



3 Dimension-Six Effective Operators 

Let us now proceed to examine in more detail the effective interaction terms 
with dimension greater than four, described schematically by the second part 
of the Lagrangian . As noted in the Introduction, we shall restrict ourselves 
to the dimension-six bosonic operators. Thus, the effective Lagrangian we are 
going to work with has a general form 

f (6) 

j? eff =j? SM + Y, J -j^of\ (10) 

i 

For the moment, let us leave aside any discussion of possible "realistic" values 
of the mass scale A and coupling constants jf^ 1 that parameterize the effects 
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of a new physics beyond SM. Here we only give a full list of the operators O- 6 ^ 
that will be utilized in our calculation. An explicit construction of such a set 
was studied systematically e.g. in [7|. Assuming that the bosonic sector of 
our effective low-energy theory contains only the particles W^, Z, 7 and H, 
and imposing the requirement of SU(2) x U(l) gauge symmetry together with 
separate invariance under charge conjugation and parity, one arrives at eleven 
independent dimension-6 monomials in the relevant field variables, namely 

Owww =Tr\W tu ,W v "W^], 

O ww = tfW^W'"'*, 

Obw = ^B^W^Q, 

O dw = Tr([D li) W vp ][D^W u "}), 

DB = -^{d ll B vp ){d»B v P), 

Obb = ^B^B^S, 
O w = {D^W^(D^), 
O b = (D^B^(D^), 

= (£» M $) t $ t $ 
0#, 2 = \d» ($t$) ($t$) , 

= I (&<S>f ■ (11) 

Some explanatory remarks concerning the notation employed in are in order 
here. For convenience, we have introduced the symbols B and W for vector field 
variables including the gauge coupling constants: 

Bpu = \ig'B^ = iffyw;, (12) 

with a a being the standard Pauli matrices and 

B^ = - amOyyZfi + cos OwA^ (13) 
W* = ±={Wl±iWZ). (14) 

The covariant derivative acting on the isospin doublet $ has the form 



B^^ + ^+ijy^ (15) 



The Higgs boson field is introduced by means of the <3? in the usual way; in the 
U -gauge this becomes 



G 



4 Contribution to H — > 77 from =5f boson 

With the basic SM results for the H — > 77 decay amplitude at hand (see 
0), a natural next goal is to calculate the corresponding leading correction 
due to the 0(1/A 2 ) terms in the effective Lagrangian Q1U[I. To this end, we 
shall consider tree-level and one-loop diagrams, in which a vertex originating in 
the appears exactly once. As we have stated in the Introduction, we are 

primarily interested in a complete evaluation of the one-particle irreducible one- 
loop bosonic graphs. Note that for this type of contributions, the transversality 
expressed by eq. © represents a non-trivial criterion of correctness of the whole 
calculation, as it is achieved only in the sum of all relevant diagrams (in contrast 
to that, a contribution of any individual reducible diagram is transverse by 
itself). 

Two remarks are in order here. First, the complete calculation would have to 
include an appropriate renormalization procedure: UV divergent parts of one- 
loop graphs, evaluated by means of dimensional regularization, should be ab- 
sorbed into a redefinition of the tree-level coefficients of the effective Lagrangian 
- in other words, they are to be cancelled by appropriate countcrtcrms, whose 
structure is determined by the considered Lagrangian HUfl . In the present pa- 
per, we do not implement such a full-fledged renormalization procedure, since 
we temporarily ignore some contributions that have been studied previously by 
other authors (for example, contributions due to dimension-six fcrmionic oper- 
ators and, in general, contributions of various reducible graphs). We prefer to 
give here just a separate treatment of the above-mentioned irreducible graphs 
(as far as we know, these have not been calculated before) and a more complete 
discussion will be presented elsewhere. For practical purposes, when display- 
ing the results obtained for the diagrams in question (see the next section), we 
remove the UV divergent parts according to the MS scheme. 

The second remark concerns the structure of the operators . Taking into 
account the £/-gauge representation of the Higgs doublet 1)16(1 . one can notice 
that e.g. the Oww ^ n d Obb contain, among other things, also parts propor- 
tional to the vector boson kinetic terms (of course, these are obtained simply by 
replacing the $ with its vacuum expectation value). In our calculation, contri- 
bution of such bilinear terms will be included in the considered Feynman graphs 
in a straightforward manner, i.e. as any other dimension-six effective coupling. 
An alternative approach would consist in incorporating such effects into the 
renormalization constants of the SM (cf. e.g. [7]). Note that bilinear interaction 
terms arise also from the other operators 

Now, let us describe in more detail the set of Feynman graphs that we are 
going to calculate. As indicated above, we first take the graphs in Fig. ^ and 
in each of them replace consecutively a SM vertex with one of those stemming 
from the Jz? b ^2 on - Similarly, one can make insertions of the 0(1/ A 2 ) terms into 
propagators within the SM graphs. In this way, one gets 16 diagrams, including 
the crossing of external legs. Next, there are several graphs with a new (i.e. not 
SM-like) topology. These are depicted in Fig. |21 (addition of the crossed diagram 
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to Fig. |3) is tacitly assumed). 




Finally, just for an illustration let us give some examples of reducible graphs. 
Two such contributions are shown in Fig. [3J As we have stated before, we will 




Figure 3: H — > 77: tadpole-type vertex corrections 



not evaluate them explicitly (though it would be relatively easy). Let us only 
make the following remark: Closed loops appearing in Fig. belong to the 
tadpole-type SM subgraphs, which are known to produce corrections to the 
Higgs vacuum expectation value v. In general, such effects can be described 
in terms of appropriate running value v(fj,) depending on the relevant renor- 
malization scale. It turns out (cf. JS]) that corrections due to running of the 
v are rather small and thus, in a complete calculation, one would not make a 
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significant error when using everywhere the basic value v(mw) = 246.22 GeV 
obtained from the Fermi constant (measured precisely in muon decay). 



5 Results and Discussion 

Let us now describe the results of our calculation. As we noted above, the 
matrix element obtained from the considered Feynman diagrams exhibits the 
desired transverse tensor structure: 

M(H - 77 ) = X(g^ - _P-) E *(p) e *(*.), (17) 

where X stands for the relevant formfactor. This can be written, in an obvious 
notation, as 

X-e 2 ff FsM + e 2 ^r + eV c ff OP (18) 
and the corresponding decay rate is then given by 

T(H -> 77) = 8 -^^^\F SM + 1*5?" + i^T ■ (19) 

Since we take into account only bosonic operators in our evaluation of the 
leading 0(1/ A 2 ) correction to the SM result, from now on we will also consider 
only the contribution of bosonic loops within SM - in other words, for the 
moment we will ignore fermions altogether. Thus, we will set 

Fsm — Fvb (20) 

henceforth (cf. (J5J). The contribution of dimension-six operators 111 l[l at the 
tree level has the following simple form: 

ptroc i ( fBB fBW fwWs 2 / 01 \ 

The non-trivial part of our calculation is represented by the closed-loop contri- 
bution. We have five basic topologies (plus two crossings); in the corresponding 
diagrams we consider insertions of effective dimension-six operators - both in 
vertices and as propagator corrections. In total, there are 23 one-loop tensor 
integrals to be calculated. Obviously, such a calculation is huge: a simple count- 
ing indicates that one has to expect tensors made of loop momentum with up 
to ten indices, schematically 



J dH I 1 " 1 ... l^ la x [denominators^ 



where six (3 x 2) Z M 's come from propagators, one Z M from SM vertices and three 
loop momenta descend from an effective vertex. Thus, the algebraic reduction 
of these expressions to the standard scalar integrals is rather long and tedious. 
Let us also emphasize that the transversality displayed in ea. (|17J) is achieved 
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only when all relevant irreducible graphs are summed; among those, only Fig.[2t 
is transverse by itself. 

The main result of this paper is the following formula for the F^ p : 

128m w n 2 F l ° s ° p = 

- 8»(/£ - 2f£ w + / A + 3g 2 (-4f£ w + rt ww ))m%m 2 w B (0,m 2 w ,m 2 w ) 

- 3(/fl B - few + fww) m H B o( m H, m H, "4) 

+ 2{m H ((/ A - 2/ A w + / A )m H - 2(/ B A - 2/ A w + ft 
+ 12s 2 (-4f£ w + ftww)W w )B Q (m 2 Hl m 2 w ,m 2 w ) 
+ [(/b ~ 2(/ A w + 7g 2 f A w ) + 3g 2 f£ ww )m 2 H - 24(g 2 f£ w + f^ w )m 2 w 

- 36(g 2 / A w + fww)r4v) Ca(m 2 H , 0, 0; m 2 w ) 
+ 4 (5 2 /bvf + fww) m w(-( m H +6m 2 v )Co(m 2 H ,m 2 H ,0;m 2 v ) 
- 2m 2 w {-m 2 H + 6m 2 w )(D {m 2 H ; m 2 w ) + 2D (m 2 H ,m 2 H ; m 2 ^)))] }. (22) 

The Bq, Cq and Dq denote the relevant Passarino-Veltman functions (see Ap- 
pendix) and we have used a shorthand notation /• = fi/A 2 . Note that UV 
divergences reside only in the Bq and can be removed by means of the MS 
prescription which amounts to the replacement 

fi = fi + til (23) 

(for the definition of the 1/e sec Appendix). Then it is not difficult to realize 
that the following relation holds 
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m 



f-(2/ B ~ if BB ~ 3f B w + Vw ~ 3/W))- (24) 



4m w 

(at least at the level of the considered irreducible graphs). 

Clearly, the result 1221) contains a lot of essentially unknown free parameters, 
so it is rather difficult to make any specific physical prediction for a correction to 
the SM decay rate. Anyway, for illustration of our general formula let us display 
at least some numerical examples. First, we discuss briefly the behaviour of the 
H — ► 77 decay width as a function of the Higgs boson mass. In Fig. 0] we have 
shown two examples of such a dependence, corresponding to different choices 
of the fi's and with A = 1 TcV. In the first option, we set all /j's equal to 1. 
The second (more or less random) choice is designed to demonstrate a possible 
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Figure 4: Bosonic contributions to T(H — > 77) for A = 1 TeV and (I): all /$ = 1 
(II): fww = 1-8, Jbb = —6.5, /bw = —3, fuw — 5-4, /www = 5.5, fs = —4 
and fw = 8. 



reduction of the T in contrast to the otherwise typical enhancement of the SM 
result. 3 

One could also wonder how much is this process dependent upon the indi- 
vidual constants /j. We could expect that the leading contribution comes from 
the tree level (cf. (|21ll ). It is of course true, but we can imagine a situation 
where the tree-level contribution is suppressed and the one-loop graphs domi- 
nate. From Fig. [3 we can see that apart from Jbb, Jbw and fww (which are 
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Figure 5: Bosonic contributions to Y{H — > 77): dependence on a given coeffi- 
cient fi with all other being set to zero, win = 120 GeV and A = 1 TeV. 



not considered in this plot because of their substantial contribution at the tree 

3 In fact, the second choice has been set so as to reduce considerably the SM result (for 
mjj = 120 GeV) by means of the tree-level contribution of dimension-six operators. 
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level), the fnw can give a rather large contribution. 

In the previous discussion we have always set A = 1 TcV, which is a generic 
estimate of the new physics scale. In Fig.[5]we have depicted the dependence of 
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Figure 6: Bosonic contributions to T(H — > 77): the dependence on A. 

the calculated width T on the A. As we would expect, effects of the new physics 
are, roughly speaking, the more important the smaller the scale A is. 

6 Conclusion 

Using the effective Lagrangian approach, we have examined the potentially in- 
teresting rare decay of the SM-like Higgs boson into two photons. We have 
considered the "decoupling scenario" , in which the onset of a new physics be- 
yond SM is supposed to be characterized by a mass scale A^> v. Within such a 
rather general framework, we have employed a full set of SU(2) x U(l) invariant 
dimension-six bosonic effective operators and evaluated, at the one-loop level, 
the leading 0(1/A 2 ) correction to the well-known SM result. Such a calculation 
seems to be missing in the earlier papers dealing with the H — > 77 process. We 
have found out that the one-loop contribution involving dimension-six bosonic 
operators can be very important if there is an accidental suppression of the tree- 
level effective interaction (i.e. when the combination Jbb — fsw + fww is close 
to zero). More generally, it is remarkable that the inclusion of dimension-six 
bosonic operators can change dramatically the bosonic SM contribution; usu- 
ally it is expected that the effects of new physics beyond SM would enhance the 
H — > 77 decay rate, but Fig. 0] shows that the effects of the operator l(TT)l could 
also reduce it significantly. 

In our calculation we have ignored completely the effect of fermions. Note 
that within SM the fermionic contribution to H — > 77 is rather small for the 
light Higgs (at least at one-loop level) in comparison with bosonic contribution 
(this follows from straightforward evaluation of (JSJ and (JBJ). Within effective 
Lagrangian approach, the effects of higher-dimensional fermionic operators were 
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studied previously in the paper ^J. Needless to say, a complete realistic cal- 
culation would have to take into account an educated guess of the values of the 
coefficients fi (based on an independent analysis of an appropriate set of other 
physical processes); in this context see e.g. |16| and references therein. Further 
work in this direction is in progress. 
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A Loop Functions 

In evaluating loop integrals we have used the so-called Passarino-Veltman re- 
duction J7|, i- e - the reduction of tensor one- loop integrals to the special scalar 
integrals which can be further expressed by means of some standard analytical 
functions. In the text we have employed the scalar integrals B , Co and D 
defined within dimensional regularization scheme as 

i^Boiplmlml) = J ( 2 ^-4 jp _ m j] [(f + Pl )2 - m 2 } ' (25) 



■ 2/~i I 2 / \2 2 2 2 _„2\ 



d d l 



Pi) 



P2) 



(26) 



iir 2 D (pl, (pi -p 2 ) 2 , (P2 ~P3) 2 ,pI,pI, (Pi ~ p 3 ) 2 ; mj, ml, ml, ml) = 
d d l 1 
(27r M ) d - 4 [I 2 - mlW + Pl ) 2 - ml][{l+p 2 y - mj}[(l + p 3 ) 2 - mj) 

The Co and Do are UV finite while Bq is UV divergent for d = 4. Defining 

1 2 

I ~ A-d 

we have 



(27) 



7b + log47r, 



1 



1-4-! 



1 



B a (p z ,m z ,m') = - - log — + 2 + Ji-4^ log- . 

One of the three-point functions used in l|22|) is given by 



(28) 



C (m 2 H , 0, 0; mw) = C (m 2 H , 0, 0; m^, m^, m^) 



•h 



Li-, 



1 - yj\ - Am 2 w /mjj 



Li; 



1 + ^\-Am 2 w /m 2 H ^ 



)], (29) 
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where hi2(x) is the standard dilogarithm defined through the Spence's integral: 

uM=-r i0g{l r t) dt. (so) 

Jo t 

Further, let us denote 

Co(m^,m^,0;m^) = C (m 2 H ,m 2 H ,0;m, 2 v ,m 2 v ,m 2 v ) , 



D (m 2 H ; mly) = D Q (m 2 H , 0, 0, 0, 0, 0; m 2 w , m 2 w , m^, m 2 ^) , 
D (m 2 H , m 2 H ; m 2 w ) = D (m 2 H , m 2 H , 0, 0, 0, 0; m 2 ^, m 2 ^, m 2 w , m^) 



A particular linear combination of these quantities appears in the last two lines 
of eq. 122(1 . The resulting expression comes out to be quite simple, namely 

2m 2 v (D (m 2 H ;m 2 v )+2D (m 2 H ,m 2 H ;m 2 v ))(m 2 H -6m 2 v )-Ca(. . .)(m^+6m^/) 
= 6(m ^ = U rctanf - " g - V (31) 



"W m H 



Note finally that the proper analytic continuation of the above functions is 
obtained by means of the prescription m 2 — > to 2 — ie wherever it is necessary. 
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